SHIP STABILITY, THEORY AND PRACTICE · VOLUME TWO · APPLIED STABILITY AND TRIM

Chapter 3 · Simpson’s Rules: Areas, Volumes, Centroids and Dynamical Stability

One arithmetic trick, four dividends: the area under the GZ curve, the area of a waterplane, the volume of the hull, and the balance points of all three.

Chapters 1 and 2 read the curve of statical stability one ordinate at a time. But the flag state, and the sea itself, care about something an ordinate cannot express: the area under the curve, which is the work the ship can absorb before she gives way. Areas under curves want integration, and the mariner’s integration is a pair of pocket sized rules published by Thomas Simpson in 1743, nearly three centuries ago. This chapter builds both rules, runs them over the GZ curve for dynamical stability and the Code’s area criteria, and then turns them loose on the hull itself: waterplane areas, the volume of displacement, and the centroids, LCF and KB, that the whole of the coming trim theory pivots on.

3.1 Why areas: the work in the curve

Heel a ship to 40° and hold her there: the moment resisting you at 40° is one ordinate’s worth of information. But to get her to 40° the sea had to work against the righting moment through every angle on the way, and work is moment times angle, summed: an area. That accumulated work is the dynamical stability, the energy a ship stores to resist external forces, and it equals the displacement times the area under the GZ curve to the angle in question. The units follow the recipe: tonnes for the displacement, metres for the lever, radians for the angle, so areas come in metre radians and dynamical stability in tonne metre radians.

10°20°30°40°50°60°70°80°0.51.0area 0° to 40°Dynamical stability= Δ × area under the curvethe work the sea must do to put her thereThe curve holds more than levers: its area is stored workan ordinate resists a push; an area resists a journeyMV Ninja, summer departure. The shaded patch is what Worked examples 3.1 to 3.3 will measure.
Figure 3.1   The shaded patch is work: the sea must spend it to put her at 40°, and she spends it back righting herself.

Degrees into radians

One radian is 57.3°, the same angle where Chapter 1 erected the GM tangent. To convert, divide degrees by 57.3: a common interval of 10° is 10 ÷ 57.3 = 0.1745 radians. Areas taken with degree intervals are meaningless; the conversion is not optional.

3.2 Simpson’s first rule: sets of 141

Slice the area into equal intervals h and measure the ordinates, the heights from the base to the curve. Simpson’s first rule fits a parabola across each pair of intervals and integrates it exactly: for three ordinates a, b, c the area is one third of h times (a + 4b + c). Chain the sets and the multipliers add where they meet: 1,4,1 becomes 1,4,2,4,1 for five ordinates and 1,4,2,4,2,4,1 for seven. The rule serves any odd number of ordinates, and the common interval must stay common: change h midstream and the rule quietly returns rubbish.

Area under curve (SR1) = ⅓ × h × (y₁ + 4y₂ + y₃)MCA formula sheet, September 2020
abchOne set of 141: a parabola on three ordinatesarea = ⅓ × h × (a + 4b + c)Stacking sets builds the multipliers1  4  1+         1  4  11  4  2  4  1+                  1  4  11  4  2  4  2  4  1five ordinates, then seven: where sets meet,the shared ordinate carries 1 + 1 = 2The rule is exact for any parabola; a fair GZ curve or waterplane is close enough to one,piece by piece, for the error to vanish in the rounding.
Figure 3.2   The first rule: a parabola on each set of three ordinates, and the famous multipliers built by stacking sets of 141.
Worked example 3.1

MV Ninja’s summer departure curve (Δ 30456 t, KG 8.09 m) computed from the booklet’s KN rows as GZ = KN − KG sin θ, gives GZ = 0 at 0°, 0.402 m at 10°, 0.866 m at 20°, 1.104 m at 30° and 1.134 m at 40°. Find the area under the curve to 40° and the dynamical stability at that angle.

HeelOrdinate GZ (m)Simpson’s multiplierProduct for area
0°010
10°0.40241.608
20°0.86621.732
30°1.10444.416
40°1.13411.134
Sum of products8.890

h = 10° = 10 ÷ 57.3 = 0.1745 radians.

Area = ⅓ × 0.1745 × 8.890 = 0.517 metre radians.

Dynamical stability = Δ × area = 30456 × 0.517 = 15746 t m r: the work, in tonne metre radians, the sea must spend to lay her over to 40°, and the energy she holds to spend back.

Laboratory 1 · The sets builder: where the multipliers come from
Each set overlaps its neighbour by one ordinate, and the shared ordinate carries the sum of the two multipliers that land on it. That single observation generates every multiplier line in the syllabus.

3.3 Simpson’s second rule: sets of 1331

The second rule fits a cubic across each set of three intervals: for four ordinates the area is three eighths of h times (a + 3b + 3c + d), and the sets stack exactly as before, 1,3,3,1 becoming 1,3,3,2,3,3,1 for seven ordinates. It serves ordinate counts of 4, 7, 10 and so on. Between them the two rules cover every count from three upwards except certain even counts, six, eight, twelve and fourteen among them, which need an extra ordinate read from the fair curve. On a fair curve the rules disagree only in the third decimal, and the examiner accepts either.

Area under curve (SR2) = ⅜ × h × (y₁ + 3y₂ + 3y₃ + y₄)MCA formula sheet, September 2020
abcdOne set of 1331: a cubic through four ordinatesarea = ⅜ × h × (a + 3b + 3c + d)Sets of 1331 stack the same way1  3  3  1+             1  3  3  11  3  3  2  3  3  1four ordinates, then seven, then ten:use it when the count is 4, 7, 10, 13...Choosing the rule: odd counts 3, 5, 7, 9 take the first rule;counts 4, 7, 10 take the second; both are welcome on a GZ curveThe two rules disagree only in the third decimal on a fair curve; the examiner accepts either.
Figure 3.3   The second rule: a cubic on each set of four ordinates. Odd counts to the first rule; 4, 7, 10 to the second.
Worked example 3.2

Using the same curve, find the area to 30°. Four ordinates are available (0°, 10°, 20°, 30°), so the second rule applies.

HeelOrdinate GZ (m)Simpson’s multiplierProduct for area
0°010
10°0.40231.206
20°0.86632.598
30°1.10411.104
Sum of products4.908

Area = ⅜ × 0.1745 × 4.908 = 0.321 metre radians.

Note the choice was forced: four ordinates cannot feed the first rule. Had ordinates at every 5° been read from the curve, seven ordinates and either rule would have served.

3.4 The three patches the Code buys

The 2008 IS Code prices stability by exactly these areas: at least 0.055 metre radians under the curve to 30°; at least 0.09 metre radians to 40° or the angle of flooding, whichever is less; and at least 0.03 metre radians between 30° and 40°, or between 30° and the angle of flooding if that is less. For MV Ninja at the summer draught the angle of flooding is 52.8°, so 40° governs. The third patch comes free once the first two are known: subtract.

10°20°30°40°50°60°70°80°0.51.0area to 30°: 0.321 m rrequired 0.055 m r30° to 40°: 0.196 m rrequired 0.03 m rwhole patch to 40°: 0.517 m rrequired 0.09 m rWhere the 2008 IS Code looks: three patches of areathe flag state buys square metre radians, and Simpson is the tillMV Ninja, summer departure: every patch clears its bar by a factor of five or more.
Figure 3.4   The Code’s three patches on MV Ninja’s departure curve, each with its bar to clear.
Worked example 3.3

Assemble the areas of Worked examples 3.1 and 3.2 into the Code’s compliance table for the summer departure condition, taking GM and the peak from Chapter 1.

Area between 30° and 40° = 0.517 − 0.321 = 0.196 m r.

Stability criterionRequirementActualCompliance
GM0.15 m2.24 mOK
GZ at 30° or more0.20 m1.104 m at 30° (maximum 1.190 m)OK
Angle of maximum GZ25°50°OK
Area to 30°0.055 m r0.321 m rOK
Area to 40°0.09 m r0.517 m rOK
Area 30° to 40°0.03 m r0.196 m rOK

Each of the five numerical criteria clears its bar by a factor of five or more, and the peak sits at twice the required angle: a deep laden bulker with a low G is a stiff, compliant ship. Chapter 15 formalises this table, adds the angle of flooding and the weather criterion’s cousins, and Chapter 16 runs it backwards into the maximum KG table.

Laboratory 2 · Raise G and watch the Code’s areas shrink
8.09 m
Summer displacement, 30456 t. The lab reads GZ ordinates at every 10° from the live curve and runs the same Simpson tables as Worked examples 3.1 to 3.3: second rule to 30°, first rule to 40°, subtraction for the band.

3.5 The same rules on steel: waterplane areas

Nothing in Simpson cares that the ordinates were righting levers. Feed the rules half breadths of a waterplane at equally spaced stations and out comes the waterplane area, doubled for the two sides. From the area, TPC follows at once, which hands us a beautiful audit: the area computed from the ordinates must reproduce the TPC the booklet prints.

0123456789101.611.912.111.92.4APFPh = 14.8 mThe ballast waterplane at 5.00 m: eleven stations, eleven half ordinateshalf breadths in red; the area is twice the Simpson integral of the half ordinatesarea = 2 × ⅓ × h × Σ(products)= 3149 m², and TPC followsFull forward, lean aft: this waterplane’s centroid sits 3.6 m forward of midships.
Figure 3.5   MV Ninja’s ballast waterplane at 5.00 m: eleven stations at h = 14.8 m, half breadths in red.
Worked example 3.4

At her ballast draught of 5.00 m, MV Ninja’s waterplane has the half ordinates below at eleven stations spaced h = 14.8 m over the 148 m between perpendiculars. Find the waterplane area and audit it against the booklet.

Station0 (AP)12345678910 (FP)
½ ordinate (m)1.67.111.911.912.112.112.112.111.911.62.4
Station012345678910
Multiplier14242424241
Product1.628.423.847.624.248.424.248.423.846.42.4

Sum of products = 319.2.

Area = 2 × ⅓ × 14.8 × 319.2 = 3149 m² (the factor of two restores the port side).

Audit: TPC = (Aw × 1.025) ÷ 100 = (3149 × 1.025) ÷ 100 = 32.28 t, and the booklet’s hydrostatic row at 5.00 m prints 32.28: the ordinates and the table describe the same waterplane, to the last figure.

3.6 Centroids: adding a lever to every product

To find where an area balances, weight each product by its distance from a chosen axis, sum the moments, and divide by the sum of the area products. Levers are conveniently counted in whole intervals from midships, forward positive; multiplying the final quotient by h converts intervals back into metres. For a waterplane the balance point is the centre of flotation, the pivot about which the ship trims, and Chapter 6 is built on it.

-5-4-3-2-10+1+2+3+4+5centroid: 3.6 m forwardThe centroid is a balance point: products of area times leverlevers counted in intervals from midships, forward positive; the fulcrum finds LCFcentroid from midships =h × Σ(moment products) ÷ Σ(area products)The same trick pointed downwards, area by area through the draughts, finds KB in Worked example 3.6.For the waterplane the balance point is the centre of flotation, the pivot of all trim in Chapter 6.
Figure 3.6   Products on a seesaw: the fulcrum that balances them is the centroid, here 3.6 m forward of midships.
Worked example 3.5

Using the ordinates of Worked example 3.4, find the position of the centre of flotation of the 5.00 m waterplane.

StationProduct for areaLeverProduct for moment
0 (AP)1.6−5−8.0
128.4−4−113.6
223.8−3−71.4
347.6−2−95.2
424.2−1−24.2
548.400
624.2+1+24.2
748.4+2+96.8
823.8+3+71.4
946.4+4+185.6
10 (FP)2.4+5+12.0
Sums319.2+77.6

Centroid from midships = h × (Σ moment products) ÷ (Σ area products) = 14.8 × 77.6 ÷ 319.2 = 3.60 m forward of midships, that is 74.0 + 3.6 = 77.6 m foap.

The booklet’s hydrostatic row prints LCF = 77.51 m: agreement within a decimetre, which is the honest price of eleven ordinates read to one decimal place. At this light draught the full bow waterplane drags the pivot well forward of midships; by the summer draught it has migrated abaft. Where that pivot sits decides how every loaded tonne changes the draughts, which is exactly Chapter 6’s business.

Laboratory 3 · The Simpson engine: your ordinates, both passes
The rule is chosen for you from the ordinate count: odd counts run sets of 141; counts of 4, 7 or 10 run sets of 1331. Edit the ordinates to any curve you please, GZ ordinates included: untick the doubling box and set h in radians for those.

3.7 Turning Simpson on his side: volumes and KB

Waterplane areas at a ladder of draughts are themselves ordinates: integrate them vertically and the answer is the volume of displacement; add draught levers and the second pass finds the height of its centroid, KB. This is precisely how the architect built the hydrostatic table, so the results must land back on it.

2400 m² at 0.0 m2939 m² at 2.4 m3132 m² at 4.8 m3319 m² at 7.2 m3442 m² at 9.6 mdraughtKB: centroid of the stackTurn Simpson on his side: waterplanes stacked into a volumeone pass sums the volume; a second, with draught levers, finds its centroid KBvolume = ⅓ × h × Σ(A products)and Δ = volume × 1.025:the stack lands on 30453 t againstthe booklet’s 30456 tEvery figure the architect prints in the hydrostatic table was born this way.
Figure 3.7   Waterplanes stacked into a volume; the centroid of the stack is B.
Worked example 3.6

The architect’s calculation sheet gives MV Ninja’s waterplane areas at five draughts from the keel to the summer marks, at h = 2.4 m. Find the volume of displacement, the salt water displacement and KB at 9.60 m, and audit all three against the booklet.

Draught (m)02.44.87.29.6
Waterplane area (m²)24002939313233193442
Draught (m)Area (m²)MultiplierProduct for volumeLever (m)Product for moment
024001240000
2.429394117562.428214.4
4.83132262644.830067.2
7.233194132767.295587.2
9.63442134429.633043.2
Sums37138186912.0

Volume = ⅓ × 2.4 × 37138 = 29710 m³; displacement = 29710 × 1.025 = 30453 t against the booklet’s 30456 t: agreement within the rounding of the areas. Note the volume itself is, to the tonne for tonne, the fresh water displacement column of the same row.

KB = (Σ moment products) ÷ (Σ volume products) = 186912.0 ÷ 37138 = 5.03 m against the booklet’s 5.041 m: within a centimetre.

The volume, the displacement and KB recovered from five areas and an eighteenth century mathematician’s rule. When the booklet, the cross curves and Simpson all tell the same story, the officer can trust the arithmetic with a ship on top of it.

Craft notes

The common interval h must stay common; if the curve ends awkwardly, read an extra ordinate from the fair curve rather than stretch an interval.

Odd ordinate counts (3, 5, 7, 9) feed the first rule; counts of 4, 7, 10 feed the second; seven suits both, and either answer earns the marks.

Levers in whole intervals keep the moment table clean; one multiplication by h at the end restores metres.

Every Simpson answer offers an audit: waterplane area against TPC, volume against displacement, moments against LCF and KB. Take the audit; it is free.

Test yourself